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Abstract 

Two integrable spin ladder systems with different types of impurities are pro- 
posed. The impurities are introduced in such a way that the integrability of the 
models is not violated. The models are solved exactly and the Bethe ansatz equa- 
tions as well as the energy eigenvalues are obtained. We show for both models that 
a phase transition between gapped and gapless spin excitations occurs at a critical 
value of the rung coupling J. In addition, the dependence of the impurities on 
this phase transition is determined explicitly. Remarkably, in one of the models a 
decreasing of the spin gap with increasing impurity strength is found. 
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The study of spin ladder systems continues to generate activity in order to understand 
the crossover from one-dimension to two- dimensions in condensed matter physics. In fact, 
with the rapid progress presently being made in nano-engineering, several compounds have 
been realized experimentally with a ladder structure (see for example ]P]>0)- In addition, 
experiments using different techniques, such as magnetic susceptibility measurements || 
or nuclear magnetic ressonance || , report on the existence of a spin gap in the spectrum 
of elementary excitations for even leg ladders. The existence of such a spin gap is crucial 
for superconductivity to occur under hole doping, as verified experimentally in some of 
these compounds and also predicted theoretically. 

Initially, most of the theoretical results concerning ladder systems were obtained from 
studies of the standard Heisenberg ladder. However, in contrast to its one-dimensional 
analogue, this model can not be solved exactly. Subsequently, in order to gain insight into 
the theory of spin ladder systems, other models with generalized interaction terms which 
guarantee integrability have been proposed |f>Hl6|l- Remarkably, such generalized models 



exhibit realistic physical properties such as the existence of a spin gap |J and the predic- 
tion of magnetization plateaus at fractional values of the total magnetization ||. However, 
although there is a considerable amount of work on integrable spin ladder systems in the 
literature, very few are concerned with the presence and influence of impurities. (One 
exception however is fllSl -) 



The role of impurities in the context of spin chains and strongly correlated electron 
systems is established to be an important aspect, particularly in low-dimensional cases. 
When undertaking studies appealing to exact solutions afforded by integrable systems, it 
is possible to incorporate impurities into the system without violating integrability. This 
can be achieved via two methods. In the first case, the impurities are given by changing 
the representation of the underlying Lie algebraic structure at some lattice sites from the 
fundamental representation to some other representation. In this context, several versions 



of the Heisenberg and t — J models, for example, have been investigated [17-20]. Another 



possibility is to introduce the impurities by way of inhomogeneities in the transfer matrix 



of the system. This was explored, for instance in |f21|-|24| . 

Here we wish to adapt these two known methods in one-dimension to incorporate 
impurities to the quasi-bidimensional case for the purpose of studying impurities in spin 
ladder systems. In particular, we will construct two integrable spin ladder models based 
on the su(4) algebra (Wang's model 0) with impurities. The effect of these impurities on 
the phase transition between the gapped and gapless spin excitations of both models will 
be investigated. We will show that in one of the cases the gap does not depend explicitly 
on the impurity, while it does in the other case. Moreover, it turns out in the latter 
instance that the spin gap decreases by increasing the impurity. 

Now, we begin by introducing the first model whose Hamiltonian reads 

N N 

JV 6 2 2 J J \ J 

H = hi ' i+1+ A(A - 8) ^' T+ A^8 hi - i+l( ^ i ~~~[fi i ~ ihLi+l+ 2 ~ ^ T '^ T ~ ^ 

i=i ^ ' " i=i ~ 



1 



where h i)i+1 = ~P i>i+1 with P i)i+1 and Q i? given by 



h iji+ i — - (1 + Of.^i+l) (1 + fi-n+i) 
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Qv = \[1 + {of of - a\a\ + o?o*)][l + (rf t* -rf^ + jf r* )]. 

Above di and t% are Pauli matrices acting on the site i (or the impurity site I) of the 
upper and lower legs of the ladder, respectively, J is the coupling constant across the 
rungs (including the impurity rung) and A is an arbitrary parameter. Throughout, iV 
is the number of rungs (equivalently, the length of the ladder) and periodic boundary 
conditions are imposed. The picture below shows the 2-leg spin 1/2 ladder in detail. 
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Notice here that besides the usual Heisenberg interactions along the legs and rungs we 
also have Heisenberg type interactions between the lattice site i and the impurity site T, as 
well as biquadratic interactions (also involving the impurity site T). In addition, there are 
three site interactions, involving the sites % and % + 1 as well as the impurity site T. Finally, 
we mention that the Hamiltonian for the ladder model based on the su(4) symmetry || 
can be recovered from eq. (0) by taking the limit A — > oo. 

The energy eigenvalues of the Hamiltonian ([]]) are given by 

E = 2 J + (1 - 2 J)N - 4 (^T " 2 ) (2) 

where \ are solutions to the Bethe ansatz equations(BAE) below. The BAE arise from 
the exact solution of the model through the nested algebraic Bethe ansatz method and 
read 



M - l_ \ -TT A; - Xj - 2l -|-r A; - flj + I 

X i + i ) IJ- + Xi-Hj-i 

l^i 3=1 J 

M 2 „. Mi , . Ms 

t-t flj - fij - Zl _ -j-r fit - Aj - I -j-r /i; - U k - I 

a i • \ Ms M 2 

Vj - A + % \ -p-r Vl - V k - 2l yr Vl- Hj-l 

v x - A - i) Al Vl - Vk + 2i Al „, -/ij + i 

The presence of the impurity can be detected explicitly through the presence of the 
parameter A in the first term of the last equation above. 

Here the ground state is given by a product of rung singlets when J > 2 and the 
energy is Eq = 2 J + (1 — 2 J)N. This is in fact the reference state used in the Bethe 
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ansatz calculations and corresponds to the case M\ = M2 = M3 = of the BAE (|3]). To 
describe an elementary excitation, we choose Mi = 1 and M2 = M3 = in the BAE (|3]) 
which gives the minimal excited state energy, E\ = —4 + 4J + (1 — 2J)N. 
The energy gap can easily be calculated and is found to be 



A = 2(J-2J. (4) 

By solving A = for J we find the critical value J c = 2, indicating the critical point 
at which the quantum phase transition from the dimerized phase to the gapless phase 
occurs. Notice that there is no effect of the impurity A on the gap. 

The integrability of this model (for J = 0) can be shown by the fact that there are 
matrices R and R* given by 
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with 



1 ,d 



2 

x 



2/x, e 

obeying the Yang-Baxter algebra 

Ru(x - y)Ri3(x)R 23 (y) = R 23 (y)Ri3(x)R 12 (x - y). 
R 12 (x - y)R* 13 (x)R; 3 (y) = i^y)^)^ - y). 

Above R(x) is the fundamental i?-matrix associated to the su(A) algebra while R*(x) is 
the solution which acts in the tensor product of the fundamental representation and its 
dual. By applying the standard procedure of the Quantum Inverse Scattering Method 
(QISM), the global Hamiltonian is obtained from 

d 



(7) 
(8) 



H=—ln(r(x,A))\ x=0 
where r(x, A) is the transfer matrix 

t(x, A) = tr (-R ,i(^) • • • R ,i( x ) R h( x ~ A ) R o,i+i(x) . . . R 0tN (x)) 
and tr denotes the trace over the auxiliary space, labelled by 0. 



(9) 



(10) 




i + 1 



N - 1 



Notice here that the impurity is incorporated into the system through the inclusion 
of the operator R*(x — A) in the transfer matrix and its effect on the spectrum can be 
detected by the presence of the extra parameter A. 
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In this way, the Hamiltonian ([!]) can be mapped to the following Hamiltonian, which 
can be derived from the R and R* matrix (eqs.(|5],|6])) obeying the Yang-Baxter algebra for 
J = 0, while for J ^ the rung interaction takes the form of a chemical potential term, 
commuting with the Hamiltonian. 



N 6 2 - - 2 - - N 

H = ) j h i<i+1 + ^ — Qi-!+ - — - h iii+ iQij- —Qi-ih iti+1 -2J y^X°° + 2JX T 00 , (11) 

i=l ^ ' i=l 

where 

3 

T _ v af3 v f3a 

a,/3=0 

Qi,i = £ xfxf. 

a,/3=0 

In the above, X± = \oti)(f3i\ are the Hubbard operators with \a>i) being the orthogonalised 
eigenstates of the local 
operator {a^.Tj). 

Let us now introduce another spin ladder model with an impurity, whose Hamiltonian 

reads 

i-l N j N j 

H = ^2 h k:k+1 + ^2 h k>k+1 + H imp + — ^2 {&k-Tk - 1) + - (ffi.Tj - 1) (12) 

fc=l k=i+l k=l 

where H imp represents the interaction of the chain with an impurity localized in the T 
position and is given by 



H imp - 4 _ A 2 



2A r -i 2A 



4- A 



(13) 

Above, hk,k+i = Pk,k+i represents the permutation operator between the lattice sites /c 
and + 1 and has the form 

= - (1 + a k .a k+1 ) (1 + ffc.Tfc+i) 

This term contains Heisenberg interactions along the legs and rung, as well as biquadratic 
interactions, whose physical importance has been pointed out in [p5|l . Here we also have 
three site interactions involving the lattice sites i, % + 1 as well as the impurity site T. 

Finally, we mention that Wang's ladder system || with length N + 1 can be recovered 
by taking the limit A = 0. 
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The energy eigenvalues of the Hamiltonian (|T2]) are given by 

2 Ml 1 J 

E = -2J+(l-2J)N + T —-iY:(wTT-o) < 14 > 



2 + A + 1 2 

where A; are solutions to the Bethe ansatz equations below 

A; - j \ N f \l - i(l - A) \ yr \l - A, - 2i -pr \l - jlj + i 

Xi + i) VAi+i(l+A)J A,-A i + 2il = lA,-^-i 

T-r /i; - /ij - 2? _ -j-r \X\ - Xj - I yr [L X ~ h> k ~ \ 

j}i M ~ Pi + 2i 7=7 M-^i + i j} = i ^i-v k + i 

M 3 . M 2 

I* vi~Vk + 2% i/, - ^ + i 

For J > (see eqs. ((L6|, [H]) below) the ground state is given by the product of rung 
singlets and the energy is —2 J + 7^ + (1 — 2J)N. This is again the reference state used 
in the Bethe ansatz calculation and corresponds to the case Mi = M 2 = M3 = for the 
BAE (p!5|). To describe an elementary excitation, we take Mi = 1 and M 2 = M 3 = in 
the BAE (|l^), which leads to the following solution for the variable A Q 



A = - 

4 24 



^{(1 + iVz)K l l* - 1(1 - zV3)(3A 2 + lG)^- 1 / 3 } (16) 



with given by 



K = 27A 3 + 12iV81A 4 + 432A 2 + 768. 



Notice that the solution of the BAE A depends explicitly on the impurity A. This relation, 
although complicated, is illustrated in Fig. 1 below. 

The energy gap can be calculated using the Bethe ansatz solution and has the form 



a = V-atttJ- (17) 

Here, in contrast to the previous case, the impurity does affect the gap. By solving 
A = for J we find the critical value J c = p 2 ^-, indicating the critical line at which 
the quantum phase transition from the dimerized phase to the gapless phase occurs. The 
phase transition line is shown in Fig. 2. 

Notice that by increasing the impurity A, the critical value J c also increases. A 
further analysis of both graphs together with the gap expression (pj]) reveals that there is 
a reduction of the gap by increasing A. This result can be easily confirmed by inspection 
of Fig. 3. 

The exact solvability of this model, as for the previous case, can be shown by the fact 
that for J = there is a .R-matrix given by equation (|5]) obeying Yang-Baxter algebra 

Rn(x - y)R 13 (x)R 23 (y) = R 23 (y)R 13 (x)R 12 (x - y). (18) 



1 strictly, the lattice length is assumed odd, such that in the limit where there is no impurity (A = 0) 
the correct gap expression discussed in § can be recovered 
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-20.0 



Figure 1: This graph shows how the solution of the BAE ( A 2 ) depends on the impurity 
A. 




Figure 2: Rung coupling J versus impurity A. This graph represents the phase diagram. 
The curve ( J c = 2/ (A 2 + 1)) divides the gapped and gapless phases. 
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Figure 3: This graph shows how the spin gap A depends on the impurity A for different 
values of the rung coupling J. 



By the standard procedure the global Hamiltonian reads 

H=^-\n(r(x,A))\ x=0 



(19) 



where r(x, A) is the transfer matrix 

t(x,A) = tr (Ro,i(x) . . . Ro ti (x)Ro ri (x - A)Ro )i+ i(x) . . . B® )N (x) (20) 



N 



N 



The shift that appears in Ro,j{x — A) represents the impurity. 

The Hamiltonian (|12|) can be mapped, as in the previous case to the Hamiltonian that 
is obtained by the standard procedure from .R-matrix (||) that obeys the Yang-Baxter 
algebra for J = 0, while for J ^ the rung interaction take the form of a chemical- 
potencial term. 

To summarize, we have introduced two spin ladder models with an impurity A . In 
the limit of vanishing impurity both models reduce to that introduced by Wang based 
on the sw(4) symmetry. The Bethe ansatz equations as well as the energy expression for 
the models were presented. We have shown that in one of the cases the impurity affects 
the gap of the system non-trivially, while in the other case there is no dependence. 
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